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Abstract 

The spectral properties of up to four interacting electrons confined within a quasi 
one-dimensional system of finite length are determined by numerical diagonalization 
including the spin degree of freedom. The ground state energy is investigated as a 
function of the electron number and of the system length. The limitations of a de- 
scription in terms of a capacitance are demonstrated. The energetically lowest lying 
excitations are physically explained as vibrational and tunneling modes. The limits of 
a dilute, Wigner-type arrangement of the electrons, and a dense, more homogeneous 
charge distribution are discussed. 



Introduction 

Interaction effects play a crucial role in the understanding of the electrical transport prop- 
erties of very small condensed matter systems at low temperatures [|IJ. Examples are 

• the Coulomb Blockade |[, where charging energies of single electrons suppress the 
current through a dissipatively shunted tunnel junction 

• Single Electron Tunneling oscillations M of the voltage across such a junction at con- 
stant current 



• resonance-like oscillations of the conductance of quantum dots, being periodic in mul- 
tiples of the elementary charge inside the dot || |6| 

An important feature in all conductance measurements on tunnel junctions and quantum 
dots is the relative isolation of the sample region from the 'external world'. In the blockade 
experiments this is achieved by a shunt impedance representing the (phase randomizing) 
influence of a coupling to the electromagnetic environment. For quantum dots weak cou- 
pling is achieved by highly resistive tunnel junctions. Here the time scale for a change of 
the electron number is large compared to the other inverse energies involved, namely the 
Fermi energy of the external wires, the charging energy and the characteristic energy of the 
dissipative heat bath. Therefore the electron number on the time scale of all relaxation 
processes approximately becomes a good quantum number. Due to its relative isolation, 
the quantum properties of the disconnected dot can be considered as one of the dominating 
factors in the single electron phenomena. The Coulomb interaction should be taken into 
account, because the charging energy is the most relevant energy scale of the problem. The 
commonly used overall description by means of a capacitance C [0, || is not completely 
obvious R [Hi and needs to be justified. This is of particular relevance, if the number of 



charge carriers is reduced to only about 10 in semiconductor samples P, ITTJ. Then both 



Coulomb interaction and kinetic (confinement) energy start to play a role. To clarify the 



interplay between these two energy contributions [|12| has been one of the main motivations 
for our quantitative study of a model of iV < 4 interacting electrons in a ID square well 
potential. We shall see, that the non-commutativity of the two corresponding energy opera- 
tors has an important influence on the ground state energy and on the excitation spectrum. 
The ground state energy in general cannot be obtained by the charging formula Q 2 /2C 
where Q is the total charge. The excitation energies have no similarity to the discrete level 
structure pertaining to the confinement. We shall present a qualitative picture, which allows 
to understand physically the spectrum in the correlated case. 

In the following Section describes our model and the calculational method. The third 
Section addresses the question of a capacitance definition in such correlated systems. One 
possibility to investigate this question is to compare the quantum mechanical ground state 
energy of N interacting, confined electrons with the charging energy of a corresponding 



capacitor. In contrast to previous work |n| [14], [15] we have been able to treat up to N = 4 



2 



electrons with high accuracy. 

The excitation spectrum is discussed in the fourth Section, for various electron densities. 
An inhomogeneous distribution is established for low charge density with essentially regularly 
spaced elementary charges as a consequence of the dominating Coulomb repulsion in that 
limit. We designate this arrangement of electrons as a Wigner molecule, in analogy to the 
Wigner electron lattice in dilute infinite systems with Coulomb interaction. Two kinds of 
elementary excitations are identified, vibrational and tunneling modes, that are characteristic 
for dilute and intermediate densities respectively. 

According to experimental situations we conclude that for the system parameters that 
characterize inversion layer based quantum dots the quantum mechanics of the electron- 
electron interaction cannot be completely discarded. In particular it should be possible to 
observe the features related to the existence of the low lying correlated excitations by suitable 
sample fabrication. 



Interacting electrons in a square well potential 
For the electron-electron interaction, we use the form 

V{x,x')<x 1 . (1) 

y/{x ~ X'f + A 2 

which behaves Coulombically at large distances. A is a measure for the width of the 
electron wave functions in transversal direction. In most of our calculations we take 
X/L — 2 • 10 -4 <C 1, where L is the system length. Then, the eigenvalues of the Hamiltonian 

H = E R ^(^H + H l ) (2) 



L \L 

depend only weakly on A . E a = e 2 /a B is the Hartree energy, a B = eh 2 /me 2 the Bohr 
radius, e the relative dielectric constant and m the electron mass. In most cases we 
present eigenvalues of the operator enclosed by round brackets. The relative importance of 
the kinetic energy in the ID square well potential 



H = J2 e n c n, CT c„ i(7 (3) 



~n< 

n,a 

( e n oc n 2 , n G IN ) decreases as compared to the Coulomb energy 



Hi — 51 Vn 4 ,nzn 2 n 1 C^ niai C^ n . ilT2 C n2(J2 C nilT1 (4) 



7ll... 714,(71, <T2 



with increasing system length L . The matrix elements V n4n3n2ni are real and do not depend 
on the electron spin a. The total spin S is therefore conserved and all eigenvalues are (2S+1) 
- fold degenerate. The symmetry relations are V n4n3n2ni = V n4n2n3ni = V n 
and V n4H3n2ni = if Ej^i = odd . For X/L < 1 

V n4n3 n 2ni ~ 87r / dq (\n(q\/L) + C) fi4(q)f 23(g) 
Jo 

where is the Fourier transform of the product .(x)<p n . (x) of the one-electron eigen- 

functions of Ho , and C = 0.577 the Euler constant. 

For the numerical diagonalization single particle states <^ na | 0> with 1 < n < M 

were chosen (usually M = 9...17 , depending on the calculation). The properly anti- 

symmetrized, non-interacting iV-particle basis ip( N ^ , including the spin degree of freedom, 

[2M\ 

is of dimensionality R = f J , 1 < v < R . In our calculations R was restricted to 

1.5- 10 4 , even when using Lanczos procedures. To avoid loops over all R 2 matrix elements of 
the Hamiltonian, we used the following economic procedure to occupy the matrix. Starting 
from the (N — 2) - particle basis, the application of two creation operators onto a certain 
V'^ V_2 ^ generates say the iV-particle state ipi with proper sign. il>( N > corresponds to 
a certain row v of the Hamiltonian matrix. Creating from the same ^" 2) a (different 
or the same !) ip^f identifies a certain column v' . The independent summation over 
all possible two-particle excitations and subsequent summation over all (N — 2) - particle 
states generates eventually all non- vanishing entries (including sums from n 4 = n\ and/or 
n 3 — n 2 1 cf. (^)) of the Hamiltonian matrix. 

Typical examples of iV-electron energy spectra are shown in Figure [TJ In presence of 
interaction, N > 1 , the density of states is much more inhomogeneous as a function of the 
energy. The lowest eigenvalues form multiplets of extremely small width when L ^> Na^ ■ 
The total number of states within each of these multiplets, including degeneracies, is 2 N . 



Ground state energies 

Figure shows the dependence of the ground state energy per particle E /N on the particle 
number iV for different L . The data are multiplied by L in order to eliminate the trivial 
L-dependence. The charging model would yield a straight line for the ground state energy 
as a function of the particle number E (N) when plotted in the same way. In very small 



4 



systems E /N deviates from a linear (N — 1) - dependence due to the discreteness of 
the spectrum of Hq . On the other hand, for systems with large L the formation of an 
inhomogeneous charge density (Wigner molecule, see below) prevents the ground state of 
few electrons energy to obey E /N oc (N — 1) . A better approximation is obtained by 
considering the Coulomb energy of N point charges at equal distances r s = L/(N — 1) . 

The importance of this charge "crystallization", which is a consequence of the charge 
quantization, for the capacity per unit length C/L can be visualized for equidistant point 
charges e in ID with a charge density distribution 

g (x) = ^8(x- Xj) , xj = ^t—tL . 

3=1 iV 1 

For a total charge Ne the capacity of such an arrangement is defined as 

C(N) := (Ne) 2 /2U 

where 



is the charging energy. For iV > 3 



L j=l ^ - J L j=2 J 



and therefore 



C{N)/L 



n 1 

1=2 J 



2(N - 1) 

In contrast to the classical capacitance of a homogeneously charged and long cylinder this 
capacitance per unit length is independent on L but explicitely dependent on the charge. 
In particular for small iV the classical concept of a capacitance obviously is inapplicable. 
Also in higher dimensionalities we expect considerable fluctuations of the capacitance as a 
function of the charge due to an inhomogeneous charge density. 

On the other hand for short systems with a more homogeneous charge distribution, 
quantum mechanical corrections to the ground state energy also do not allow the use of the 
capacitance formula. In Figure |3| E /N is shown as a function of r s . The deviations from 
the capacity like 1/L behavior occurring below r s <, 50cib cannot be attributed to a simple 
additive influence of the kinetic part of the Hamiltonian. The ground state energy, of the 
non-interacting system has already been subtracted in Figure [| 



Low lying excitations 

For L ^> Nob the spectrum of the low lying excitations can be understood using the 
picture of a Wigner molecule. The one particle density shows N approximately equidistant 
peaks |TB |. One type of excitation in such an arrangement is of a phonon kind due to the 
Coulomb forces between the charges. Similar to the one particle density these excitations 
are insensitive to the total spin and the symmetry properties of the wave function. 

To estimate the asymptotic behaviour of typical phonon frequencies Q as a function of 
the electron distance r s , we assume ^-functions or Gaussians for in the charge density of 
each peak. The Gaussians can be considered to emerge from the harmonic oscillator ground 
state wave functions due to the linearized electrostatic potential for one certain electron, 
leaving the other electrons fixed. The result of this crude estimate is a r~ 7 decrease of Q 
with 7 = 3/2 (5-function density) or 7 = 1 (Gaussian density). Figure |] shows VtL/a-Q as 
function of r s for different N , Q is the distance between the lowest two multiplets in our 
spectra. The behavior at large r s is consistent with 7 >, 1 . This indicates that the charge 
density distribution is more localized than a Gaussian. For r s < lOOae the strong deviations 
from the asymptotic behavior signalize the breakdown of the Wigner molecule. 

The multiplets consist of a series of energy levels which, compared to Q , are extremely 
close to each other. Each level is (25 + 1) - fold degenerate since all eigenstates of H are 
simultaneously eigenstates to S 2 . For N > 3 the states are in general not products of 



a spatial part and a spin part \ Y7j . The Pauli principle only requires their transformation 
according to the totally antisymmetric irreducible representation of the symmetric group 
Sn ■ 

In the following a qualitative understanding for the fine structure spectrum within one 
multiplet is developed, starting from the notion of the Wigner molecule. Quantitative results 
will be published elsewhere JTB ]. The total potential, including the interaction ([TJ), has N\ 
equivalent minima in the configuration space which is a iV-dimensional hypercube L N . 
There are no further equivalent minima, because Sn is the only symmetry group of the 
problem. For sufficiently low electron densities (the corresponding r s will be quantified by 
numerical results below), the probability amplitude of the eigenfunctions is well localized 
around the vicinity of the potential minima. The barrier of the Coulomb interaction in the 
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ID case separates adjacent minima. 

Each of the probability amplitude peaks may be taken as one of the basis functions of a 
finite dimensional "pocket state basis" |TS| {|j>} , with 1 < j < N\ . Matrix elements of 
the Hamiltonian in this basis <j \ H \j'> = Hjji mutually connect different permutations of 
iV-electron states and behave like tunneling integrals. The smallness of the Hjji is crucially 
required to approximate the true interacting eigenf unctions of H by eigenvectors of the 
Hamiltonian in this truncated basis. They are given by certain linear combinations of the 
{|j>} . Correspondingly, the true energy eigenvalues are approximated by eigenvalues of 
the pocket state Hamiltonian matrix. Since the amplitudes of the \j> decrease roughly 
exponentially with the distance from its center, the Hjji have also this property and thus 
the pocket state approach improves with increasing r s . It is important to note, that this 
description is not based on any single particle basis set as a starting point. It becomes in- 
creasingly reliable with increasing influence of the interaction. However, the diagonalization 
of the Hamiltonian matrix in this basis cannot be achieved by Fourier transformation, since 
Sn is not abelian. 

As a further approximation, we neglect all Hjji except those which correspond to nearest 
neighbours in the configuration space L N . For sufficiently large r s this is justified by the 
exponentially fast decay of the pocket state far away from its center. Due to the symmetry 
of the problem, all of the remaining tunneling integrals Hjji = t must be equal and are 
the only nontrivial entries in the Hamiltonian matrix. Therefore the differences between the 
eigenvalues scale with the common factor t . 

In Figure |5] the energy difference A between the two lowest numerically obtained eigenval- 
ues is plotted versus the system length L . The assumption of an asymptotically exponential 
decay of the pocket states suggests that A oc exp(— L/L&) with La ~ 1.5a B from Figure [5| 
This establishes a length scale beyond which the spectrum of non-interacting electrons is 
changed into one showing the multiplet structure that is characteristic for the influence of 
the Coulomb interaction. The tunneling energies A depend on A , because the height of 
the barriers between the potential minima are proportional to A~ x (cf. (p])). 

Conclusion 

We have calculated numerically the energy spectra of up to N = 4 electrons confined in a 
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quasi-one dimensional square well potential of finite length. The discussion in terms of the 
pocket state basis suggests that our classification of the energy eigenvalues should remain 
valid also in 2D or 3D systems, if the width of the system does not exceed the width of one 
pocket state wave function. For iV = 2 we can reproduce the fine structure in the lowest 
multiplet calculated by Bryant |L3] with its sequence 5 = 0,1,0 of total spins and with 
equal level distances for sufficiently large size of a rectangular two dimensional quantum dot, 
which is 10 times longer than wide. 

We have demonstrated that the ground state energies E (N) deviate from the N 2 - 
behavior assumed in the charging model because of the formation of a Wigner molecule like 
structure at sufficiently low electron densities (quantization of the charge) and the quantum 
mechanical influence of the kinetic energy (non-commutativity of H and Hi ). Only 
in sufficiently large systems and at sufficiently high electron densities, a capacitance-like 
behavior can be obtained. 

We have obtained 3 different regimes for the electron densities to characterize the excita- 
tion spectra. The Wigner molecule is found to be fully established for densities l/r s below 
lO^ctg 1 (Figure f|, see also |16|])- Nevertheless, the description of the interacting spectrum 
in terms of the pocket state picture does already hold at much larger electron densities. Only 
for L^L^ m 1.5cib the confinement energy starts to dominate the Coulomb energy and 
the spectrum approaches the non- interacting limit (Figure |5|). Neither the ground state 
energy nor the level spectrum is given by a sum of kinetic and potential energy eigenvalues 
separately ! 

Experiments are frequently performed on AlGaAs/GaAs - based heterostructures which 
rather correspond to a 2D situation. It is not obvious in how far our ID classification for the 
length- and energy scales of few Coulombically interacting electrons can be applied to that 
case. If we assume that at least the qualitative aspects of our classification into different 
regimes for the electron density can be used, the intermediate regime should apply in most 
circumstances. Given the geometry and the electron numbers in typical quantum dots || 
(area of the dot ~ 10 5 nm 2 , number of electrons ~ 10 2 effective mass ~ 0.07m , dielectric 
constant w 10) a mean distance of r s « 3a# can be estimated. For this relatively high 
electron density and number the ground state energy, which is in first approximation the 
relevant quantity that enters a dc-transport experiment, can roughly be estimated by using 
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the charging model. However, the excitation energies are qualitatively different from the 
ones expected within the non-interacting picture. They are importantly characterized by 
the fine structure level spectrum and the total spin. In experimental situations as they have 
been realized recently by Meurer et. al. JTl|] with only a few electrons per quantum dot, the 
charging model even cannot be expected to yield correct results for the ground state energy. 

It should be possible to observe the excitation spectrum obtained in this paper in opti- 
cal measurements |2^, |TTJ (if the potential is not strictly harmonic |15||) and in non-linear 
transport measurements [O at low temperatures. At least for quasi ID geometries, where 
the spatial width of the dot region is of comparable size or smaller than the width of the 
pocket state wave function, we predict a close relation between the energies of low lying 
excitations and the number of electrons in the dot. Only a finite number of tunneling type 
excitations should exist for fixed electron numbers. The energy scale for a 100 nm long ID 
structure on AlGaAs/GaAs basis can be estimated from Figure to be A « 0.1 — 0.5 meV 
for N = 2 . . .4 . However, this value still depends on the effective width of the structure. 

For the non-linear transport measurements two remarks should be made, i) the excita- 
tion spectra of N- and (N + 1) - electron states are completely different, even the typical 
splittings may differ by almost an order of magnitude. They are not equidistant and inde- 
pendent on N . ii) the total spin of the electrons inside the dot region can be changed only 
by ±1/2 when an electron is added or removed at the finite conductance situation. This 
reduces the number of transport channels available. For example the conductance channel 
connected with the energy difference ^5=3/2 — Eg^ 2 ' ( /i^N = 2) in the nomenclature of 



Ref. [12], S is the total spin) should not appear. 
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Figure Captions 



Figure 1: Typical spectra of model (^) for various N and L = 9.45ae • For 
N > 2 the low lying eigenvalues form groups of (fine structure) multiplets, 
the total number of states per multiplet being equal to the dimensionality of 
the spin Hilbert space 2 N . For clarity the lowest multiplets are magnified 
indicating the total spin of each level. The ground state energy is subtracted 
respectively. 

Figure 2: Ground state energies per particle E /N multiplied by L/a B versus 
the particle number N for L = 6.61ae (□) L = 16.1ob (°) L = 94.5ae (A) 
L = 944. 8ae (+)• ( x ) denote the energy of iV fixed point charges equally 
spaced at distances Lj (N — 1) . The quantum mechanical ground state energies 
(slowly) approach these values as L — >• oo . 

Figure 3: Ground state energy per particle multiplied by L/aB as a function 
of the electron distance r s := Lf(N — 1) ; the corresponding non-interacting 
ground state energy is subtracted. Pronounced deviations from the Coulombic 
l/L behavior occur below r s ^SOae • These deviations cannot be attributed 
to a simple additive influence of the kinetic part of the Hamiltonian. 
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Figure 4: Energy difference Q between the two lowest multiplets, multiplied 
by L/a-B versus the mean particle distance r s for N = 2 and N = 3. For 
r s >, lOOae the asymptotic behaviour is recovered. 

Figure 5: Logarithm of the energy difference A between the ground state 
and the first excited state within the lowest multiplet versus the system length 
for N = 2, M = 11 (□), N = 3, M = 13 (o), and N = A, M = 10 (A). From 
the slope we estimate La ~ 1.5ob. 
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